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Acoustic cloaking in 2D and 3D using finite mass 
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Fundamental features of rotationally symmetric acoustic cloaks with anisotropic inertia are de- 
rived. Two universal relations are found to connect the radial and transverse phase speeds and 
the bulk modulus in the cloak. Perfect cloaking occurs only if the radial component of the density 
becomes infinite at the cloak inner boundary, requiring an infinitely massive cloak. A practical cloak 
of finite mass is defined in terms of its effective visible radius, which vanishes for perfect cloaking. 
Significant cloaking is obtained when the effective visible radius is subwavelength, reducing the total 
scattering cross section, and may be achieved even as the interior radius of the cloak is large relative 
to the wavelength. Both 2D vs. 3D effects are compared as we illustrate how the spatial dependence 
of the cloaking parameters effect the total cross section. 

PACS numbers: 43.20.+g,42.25.Fx 



The fundamental observations of Pendry et al. [l[ and 
of Leonhardt [2J that the electromagnetic equations re- 
main invariant under spatial transformations has gener- 
ated significant interest in the possibility of passive acous- 
tic cloaking. The idea is to map a region around a single 
point in such a way that the mapping is one-to-one ev- 
erywhere except at the point, which is mapped into the 
cloak inner boundary. 

Consequences of the transformation method for EM 
cloaking have been studied extensively. Thus, Rahm et 
al. [|| demonstrate square cloaks; Chen and Chan [i| pro- 
vide an interesting mapping that rotates the EM field in- 
side the cloaking region; Chen et al. @ examine in detail 
the interaction with a spherical cloak characterized by 
the linear radial transformation Zhang et al. @ con- 
sider cloaking in inhomogeneous environments; Huang et 
al. 0] describe a method to realize cylindrical cloaking 
using a structure with "normal" EM layers; Kwon and 
Werner Q extend the linear transformation to an eccen- 
tric elliptic annular geometry lacking rotational symme- 
try; a quadratic cloak was proposed by and investigated 
by Yan et al. @; Cai et al. [l(| use a similar quadratic 
spatial transformation to obtain a smooth transition in 
the moduli at the outer interface; and Leonhardt and 
Philbin pd| provide a general theory of the transforma- 
tion method in the context of EM waves. 

An acoustic fluid is described by a bulk modulus K and 
density p, but cloaking cannot occur if the bulk modulus 
and density are scalar quantities. It is possible to ob- 
tain acoustical cloaks by assuming a two parameter den- 
sity tensor [l2|, EH, E3 ■ A tensorial density is not ruled 
out on basic principles fl5j and may in fact be realized 
through so-called metamaterials. The general context for 
anisotropic inertia is the Willis equations of elastodynam- 
ics [l|| which Milton et al. [HI showed are the natural 
counterparts to the EM equations that remain invariant 



under spatial transformation. Acoustic cloaking has been 
demonstrated, theoretically at least, in both 2D and 3D: 
a spherically symmetric cloak was discussed by Chen and 
Chan (l3j and by Cummer et al. [12], while Cummer 
and Schurig [3] described a 2D cylindrically symmetric 
acoustic cloak. These papers use a linear transformation 
based on prior EM results in 2D [llj. While the focus 
here is on passive cloaking, we note that Miller [lj| de- 
scribes possible active acoustic cloaking strategics. 

The 2- and 3D cloaks of [H, EH, E3 are shown here 
to be special cases of a more general class of acoustic 
cloaks. We will show that the linear transformation pro- 
posed in [l[ and examined by [H, [l2j| is just a special case 
of a general class of transformations possible in acoustic 
cloaking. The arbitrary nature of the spatial transforma- 
tion function has been noted flfjj but not utilized much 
as yet, except for quadratic cloaks [§, E(|. In this paper 
we show that the cloak density tensor depends upon the 
arbitrary spatial mapping function. Issues that have not 
been previously considered, such as the total amount of 
mass required for cloaking, are addressed. It turns out 
that perfect cloaking occurs only if infinite mass is avail- 
able, regardless of the spatial mapping function. We take 
the point of view that achieving infinite mass is imprac- 
tical, but that effective cloaking with finite mass is still 
possible if the cloak parameters are appropriately chosen. 
We begin with a derivation of the general form of the 
density tensor and bulk modulus necessary for acoustic 
cloaking with rotational symmetry in both 2D and 3D. 

The equations governing small amplitude disturbances 
in an acoustic mctamatcrial are 



pv = — gradp, 



-K div v. 



(1) 



Here, p(x, t) and v(x, t) are the pressure and particle ve- 
locity field variables, and if(x) > is the bulk modulus. 
The density, p(x), is a symmetric second order tensor. 
Traditional acoustics corresponds of course to the diago- 
nal form p = pi. Equations |l| imply that the pressure 
satisfies the generalized acoustic wave equation 
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K div (p 1 gradp) — p = 0. 



(2) 
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We examine the 2D and 3D configurations in parallel. 

Consider spatially inhomogeneous anisotropic materi- 
als that are locally transversely isotropic, that is, charac- 
terized by an axis of symmetry and an orthogonal hyper- 
plane of isotropy. The axis of symmetry is in the radial 
direction, such that the density matrix has the form 



p = p r (r)x Cg) x + p±(r)(I - x <X> x, 



(3) 



where r — |x| and x = x/r. The bulk modulus, being 
a scalar, depends only on radius, K(r). The governing 
equation for p is 



K(r) d 



' r d 1 dp~\ | K(r) 
^p r {r) dr 



■A ± p-p = 0, (4) 



where d is the dimension (2 or 3) and Aj_p is the 
Bcltrami-Michel operator. 

The central idea is to contract or shrink the radial 
coordinate r — > f(r), with the pressure parameterized 
in terms of the contracted radial coordinate as p(x, t) = 
P(/(r), x, t) . The scalar governing equation ([4]) becomes 



K(r) 8 (r*~\ 



f dP\ K{r) 



df\Pr(r) df 
where /' = df/dr 



r 2 p± 



A ± P-P = 0, (5) 



equation is V^P — P 



At the same time, the scalar wave 
= in the coordinates (/, x) where 
W'j is the Laplacian in these un-contracted coordinates. 
The pressure simultaneously satisfies this uniform wave 
equation in (/(r),x) and eq. ([5]) in (r, x) if and only if 
the following three conditions are met: 



r.d-1 



„d-3 



r.d-1 



f 



f 

d-l J ' 



K 



(6) 



The exterior of r = b is a uniform and isotropic acoustic 
medium with K = 1, p — I. Continuity at r = b requires 
matching of the pressure and the normal (radial) velocity, 
and hence the radial acceleration. The latter follows from 
eq. ((T|) as v r = —p~ 1 dp/dr, or v r = —(f/r)f d ~ 1 dP/df. 
Hence, all that is required for continuity at r = b is that 
/ be continuous across the boundary, which is accom- 
plished by requiring f(b) = b. Surprisingly, the two con- 
ditions at r — b reduce to this single requirement. 

Equations ((6|) provide relations for the anisotropic 
acoustic material properties for a given contraction /(r). 
It is useful to express them in terms of the radial and az- 
imuthal phase speeds: c r = yf Kj p r and c±_ = yj Kj p± . 
The mass density tensor can then be expressed p = 
Equations j6|) imply the identity 



K 



cT 2 !, 



K 



PrPL 



d-l 



(7) 



independent of the choice of contraction function f(r). 
The quantity Kp r is the square of the radial acoustic 
impedance, z r = Kp r . Equation implies that the 
identity z r = c^T 1 is required for cloaking. The three 



equations ([6]) can be replaced by along with relations 
for the phase speeds in terms of /: 



=r = 1//', c x = r/f. 



(8) 



Note that /' is required to be positive. The original quan- 
tities can be expressed in terms of the phase speeds as 



c 



p± = C r C i 



d-3 



K = c 



„d-l 



r 



(9) 



One could, for instance, eliminate / as the fundamental 
variable defining the cloak in favor of c± (r) , from which 
all other quantities can be determined from the differ- 
ential equation relating the speeds: (r/cj_) = l/c r . We 
assume the cloak occupies H = {Q<a<r<b} with uni- 
form acoustical properties K = 1, p = I in the exterior. 
Then cj_(r) follows by integration, 



1 



b d 



— , a < r < b. 



(10) 



Can the cloak density be isotropic? This will occur if 
c r = cj_, which requires that /' = f/r, implying / = jr 
with 7 constant. Impedance matching at the cloak outer 
boundary r = b requires f(b) = b. This can only be 
satisfied in the trivial case of 7 = 1, implying f = r. 
While this eliminates the possibility of isotropic cloaks, 
it is perhaps a satisfying outcome. 

Perfect cloaking requires that / vanish for some r = 
ao > 0. It is clear that c± and z r blow up as r J. ao- 
Hence the product Kp r blows up, but do K and p r each 
become unbounded? Consider / oc (r — ao) a near ao 
for a constant and non-negative. No value of a > will 
keep the radial density p r bounded, but the unique choice 
a = 1/d ensures that the bulk modulus K(ao) remains 
finite. We will see below that this power law is not a 



R=f(r) 



R=r 
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FIG. 1: The cloak occupies a < r < b with uniform 
properties outside and the " 'hidden" ' object inside. The 
cloak properties depend upon a monotone contraction 
mapping /(r) within the cloak, such that f(b) = b, 
f(a) < a. It is useful to think of an invisibility 
boundary r = ao inside the cloak, where /(ao) = 0. The 
pressure satisfies a uniform wave equation in the 
mapped radial variable R(r) . 
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(a) 2D, a = 1, a ta t = 1.18 (b) 3D, a = 1, <r toi = 0.49 




(c) 2D, a = 4, <r tot = 0.06 (d) 3D, a = 4, tr toi = 10~ 4 

FIG. 2: A plane wave is incident from the left with 
frequency k = 10 on the cloak defined by eq. (fT3"| with 
(&, a, ao) — (1,0.5,0.35). The exponent is a = 1 in (a) 
and (b) implying a virtual inner radius of /(a) = 0.23. 
For (c) and (d), a = 4 giving /(a) = 0.003. 



TABLE I: Behavior of quantities near the vanishing 
point r = ao for the scaling / oc x a as x = r — a I 0. 
The total radial mass m r is defined in eq. (jlljl . 
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practical choice for the cloak mapping function. Note 
that the azimuthal density has a finite limit in 2D for 
power law decay / oc (r — ao) Q , while p± remains finite in 
3D iff a < 1 , otherwise it blows up. Similarly, the radial 
phase speed scales (r — ao) 1 a , which remains 



finite for a < 1, blowing up otherwise. These results are 
summarized in Table I. 

We introduce a nondimensional measure of the total 
amount of cloaking mass present in the radial component 
of the inertia tensor: 

J n dV b d -a d J v ; 

a 

Explicitly, 

f wh? [ b2 ln /( & ) - « 2 ln /(«) - 2 /' dr r In /(r)] , 
m r = < 

(12) 

in 2D and 3D, respectively These forms indicate not 
only that m r — > oo as /(a) — > 0, but also the form of the 
blow-up. To leading order, m r = b ?°i a -} In j^y + ■ ■ ■ i n 

o 3 

2D and m r = b :i„ a a j^y + • • • in 3D. The blow-up of m r 
occurs no matter how / tends to zero; that is, the infinite 
mass is an unavoidable singularity. 

The idea behind the contraction function /(r) is to 
map the entire interior of r < b into an annular subset: 
Q. This in turn implies that perfect cloaking requires 
infinite mass m r . We take the point of view that m r must 
remain finite, which can be achieved by making / > in 
Q. To be specific, consider the analytic continuation of 
f(r) into the "cloak" {r : < r < a}, where wc allow f(r) 
to vanish. Consider the power law form for the mapping 
function, f(r) = f a (r) where 

/ (a) M^(P) a . (13) 

for < ao < a. This yields finite mass m r , but it blows 
up as a | ao for every a > 0. 

It follows from (fT2|) that the dependence of the total 
mass on a is linear in 2D: rrS^ = araf'. In other words, 
in 2D the mass required is proportional to the log of 
the effective vanishing radius. For 3D, the leading order 
term in eq. (|12p implies w (3~ 1 (/3mf 1 '') Q , where 

(3 = H (|y — 1) . The total mass grows exponentially with 
a, all other parameters being fixed. However, we can still 
achieve effective cloaking with finite m r . The idea is to 
make /(a) subwavelength. This does not require that a 
be subwavelength. 

We assume time harmonic motion, with the factor 
e -ikt unc j crs tood but omitted. Consider a cloak of finite 
mass with zero pressure on the interior surface r = a. 
The total response for plane wave incidence is 
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p = p e 



ikR(r) cos 6 



^"(2-^ )J„(fcfl(a)) ggggj cosnfl, 2Z>, 



»"(2n + l)j w (fefl(o)) ffiggg ^(cos g), 3D, 



f(r), a<r<b, 



R(r) 



(14) 



r, 



b < r < oo. 



The total scattering cross-section, and hence the to- 
tal energy scattered, depends on the forward scatter- 
ing amplitude through the optical theorem: a tot = 
47r/c _1 Im<?(e 2 ), where g defines the leading order far- 
field p = p e lkz +p g(Z)r-( d -V/ 2 {i2Tr/k)( 3 - d ^ 2 e ikr . The 
scattering cross-section is easily computed from (1141) , and 
is dominated in the small kR(a) limit by the n = term, 
with leading order approximations 

r^|mfc/(a)|- 2 + ..., 2D, 
otot = < (15) 
[4Trf 2 (a) + ..., 3D. 

The faster decay of <Jtot explains the greater efficacy of 
the finite mass cloak in 3D. Considered as a function of 
the power a, in 3D the mass m r increases exponentially 
while the cross-section a tot decreases exponentially In 
2D the rates of increase and decrease of m r and Otot are 
algebraic (linear) or less. The comparison suggests that 
cylindrical cloaks require that the visibility radius /(a) 
be very small. Figure 2 illustrates clearly the disparity 
between chylindrical and spherical cloaking. Thus, for 
f(a)/b = 0.003 the 3D cross-section is negliblc (Fig. 2.d) 



but for 2D the cross-section is two orders of magnitude 
larger (Fig. 2.c). Note that Ruan et al. [3 found that 
the perfect cylindrical EM cloak is sensitive to perturba- 
tion. This sensitivity is evident from the present analysis 
through the dependence on the length a — ao > which 
measures the departure from perfect cloaking (a = ao). 

In conclusion, we have shown how to generate acoustic 
cloaks in 2D and 3D using a general radial transrforma- 
tion r — ► f(r). Regardless of the form of /, two universal 
equations relate the cloak parameters, e.g. K — c r c^ _1 
and (r/cj_) = l/c r . Perfect cloaking requires not only 
that the density p r — > oo at the inner radius, but the 
associated total mass m r is also infinite. This can be 
avoided in practice since the scattering cross-section a tot 
can be made as small as desired by selecting the value 
of f(a) <C a so that kf(a) <C 1 even as ka =0(1). We 
showed that the form of the cloak mapping function f(r) 
is important and the effect of the power a is markedly 
different in 2D vs. 3D. The feasibility of cloaking de- 
pends on the development of real metamaterials with the 
desired unusu al phys ical properties. This is an active 
research area [2(J |2l|, [22j with tremendous potential for 
future applications. 
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